MASSACHUSKITSINSTrrUTROFTKCHNO'LOGY 

AR'I 1F1CIAI INTELLIGENCE LABORATORY 

A. I. Memo No. 612 January 1981, revised March 1983 

THE CURVE OF LEAST ENERGY 
B. K. P. Horn 



^""N 



Abstract: Here we search for the curve which has the smallest integral of flic square of curvature, while 
passing through two given points with given orientation. This is the true shape of a spline used in lofting. In 
computer-aided design, curves have been sought which maximize "smoothness". The curve discussed here is 
the one arising in this, way from a commonly used measure of smoothness. The human visual system may use 
such a curve when it constructs a subjective contour. 
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Introdiiction 

The curve which passes through two specified points with specified orientation while minimizing 

8 = §k z ds, 

where k is the curvature and s the arc-distance, has a number of interesting applications. 

In a thin beam, curvature at a point is proportional to the bending moment [pg. 80, 1]. The total elastic energy 
stored in a thin beam is therefore proportional to the integral of the square of the curvature [pg. 163, 11 The 
shape taken on by a thin beam is the one which minimizes its internal strain energy. This is why we call the 
curve sought here the minimum energy curve. A thin metal or wooden strip used by a drafts(wo)man to 
smoothly connect a number of points is called a spline [pg. 156, 2]. Such splines are used in creating lofted 
surfaces from plane parallel cross-sections of ship hulls and aircraft fuselages [pg. 228, 2]. The shape of a 
spline constrained to pass through two specified points with specified orientation is what we are after here. 

In computer graphics and computer aided design there is a search for curves which are particularly "smooth" 
[pg. 156, 2; pg. 49, 3; pg. 119, 4; pg. 309, 5; pg. 43, 6]. One measure of smoothness is the inverse of the 
integral given above. Typically, cubic polynomial approximations are used instead of the optimal curve 
[pg. 162, 2; pg. 66, 3; pg. 129, 4; pg. 315, 5]. (Unfortunately, these approximations are called splines too). 

It has been suggested that the human visual system uses a curve of low energy when completing a contour. 
Ullman [pg. 1, 7] proposes that a subjective contour consists of two circular arcs tangent at their point of 
contact. Out of the one-parameter family of solutions of this form he picks the one which minimizes the 
integral of the square of curvature [pg. 2, 7]. He notes that the curve so constructed may have near-minimal 
energy. (This does not mean that it necessarily lies close to the curve of minimum energy, as we shall see.) 
Brady ei al used cubic polynomial approximations instead [8]. 



Preview 

We will first consider a special case. Here the curve must pass through the points (-1,0) and ( + 1,0) in the 
xyphnc with vertical orientation at both points. We first determine the optimum curve as the limit of a series 
of approximations. This approach helps suggest algorithms for computing approximations to the ideal curve. 
Later we solve the variational problem directly. We then discuss how this curve can be translated, rotated, 
and scaled to produce a four parameter family of curves which contains the general solution. 



A Semicircle 

One curve which connects the two points (-1,0) and ( + 1,0), and has the desired orientation at these points, is 
a semicircle of radius one with the center at the origin (see Figure 1). The curvature equals one at all points 
and so the relevant integral for the section in die right hand quadrant becomes 

8= f* /2 l& = */2. 

Jo 

The arc length happens to have the same value, 






r\ 
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whilc the maximum height of the curve above the x-axis is 

% = 1. 
Can we do better, that is, find a curve with a smaller value of 8? 

Two Circular Arcs 

We try a combination of two circular arcs for the portion of the curve in the first quadrant, as shown in 
Figure 2. (The other half of the curve is obtained by reflection about the >'-axis.) Let the first arc have radius 
R and angular extend 0, while the remaining portion has radius rand angular extend (w/2-0). Note that the 
parameters t\ /?, and 6 are not independent, since one can obtain from the diagram 

(R-r) cos 6 = (rt-1). 

The arc length of the right-hand portion of the curve becomes 

y = R + (tt/2 -0) r, 

while die energy is 6 = 1 + ^ /2 "^ . 

R r 

Minimizing g, subject to the given constraint, using the method of Lagrangian multipliers, leads to the set of 
equations 

^ + ^(l-sec*) = 0, 

(/?-/-) cos -(R-l) = 0. 
The second of these can be simplified into 

r = R (n/2 d) tan 6, 
which, when applied to the first equation, yields 

6 (tt/2 -0) (1 + sec 6) = 1. 

Solving this numerically we obtain 8 - 0.412868765... 

We can also show that R - rp. .. , \ .. .. . . - , , 

[(1 -cos 0) + (n/2 -0) sin 0] ' 

so that R = 1.8227161... 

and r = 0.92452847... 
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Finally, since here % = r + (R - r) sin B 

wc see that DG = 1.2849161... 

and 6 = 1.4789649... ~ 0.94153834*(tt/2), 

wh ile *= 1.8230795... = 1.1606084*(w/2). 

The energy in this curve is only 94.15% of that in the semicircle, so we can do better. 

The Best Ellipse 

The two-arc solution suggests that the optimum curve is elongated and has radius of curvature smaller than 
one at its peak, and larger than one near the jo-axis. Is it an ellipse? The equation of the ellipse [pg. 411, 13; 
pg. 72. 14] shown in Figure 3 is 

(x/a) 2 + {y/b) z = 1, 
or in parametric form, x = a cos /, y=b sin /. 

The eccentricity e is defined by the equation 

e 2 = 1 - (a/bf. 
The arc length can be found as follows [pg. 26, 9], 

,= /*= JT' 2 V * +fdt 

where k= d mdy = ~df 

x 2 + y 2 = a 2 sin 2 / + b z cos 2 ?, 



jc 2 + f = b 2 [1 - e 2 sin 2 /]. 



Now, r + y 

or, 

So 3=b f sfl-e 2 sin 2 / A = 6 £(e)> 

where E(e) is the complete elliptic integral of the second kind 1 [pg. 16, 9; pg. 833, 10; pg. 904, 11; pg. 589, 12]. 
(The elliptic integrals got their name from the fact that they first appeared in the mensuration of the ellipse.) 
Finally, for a-\, 

X = E(e) I -y/ \-e z . 



1. Note that E denotes two things in this paper. When it has one arguments here, it signifies the complete 
elliptic integral of the second kind. 
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fhc curvature can be found as follows [pg. 553, 13; pg. 22, 14]: 



K , yi-yx ab 

[^+;' 2 ] 3/2 "b 2 sin 2 /+6 2 cos 2 /] 3/2 - 



Now 6 = JV ds = J 77 ' 2 k 2 [i 2 + /] 1/2 <*, 

2 cn/l 



bO b = 73 I 77 2 • 2 i5/2 "'• 

6 d Jo [l-6 z sin z /] 5/z 

r Fhe definite integral can be shown [using pg. 165, 11] to be equal to 

W l + \)E{e)-e l K(e)\IQe*l 

where e 2 + e 2 = 1, 

and #(<?) is the complete elliptic integral of the first kind [pg. 16, 9; pg. 834, 10; pg. 904, 11; pg. 589, 12]. So 
finally when a= 1, 

8 = [2 (2 -e 2 ) E(e) - (1 -e 2 ) K(e)] / [3 ^ l-e 2 \ 

When we set e=Q we obtain 8 = tt/2, as we should, since the curve in this case is just a semicircle. 

In order to find the best possible ellipse we need to differentiate the expression above with respect to e. Here 
we need the following derivatives [pg. 21, 9; pg. 907, 11]: 



dE{k) 

dk ~ 


E(k) 


-K(k) 
k 


dK(k) 


E(k) 


- k' 2 K(k) 


dk ~ 




kk' 2 



where k 2 + k' 2 = 1. 

The eccentricity of the optimum ellipse satisfies 

E(e) [4e A - Se 2 + 3] = K(e) [2e A - 5e 2 + 3]. 
Solving this equation numerically leads to 

e = 0.6530018... 
with %-b/a-ll sfT^J = 1.3203823... 

8 = 1.4674751... 22 0.93422368 * (w/2), 



y = 1.8311202... « 1.1657273 * (tt/2). 
The maximum and minimum radii of curvature are 

r = b 2 /a = 1.7434096... and r. = <///> = 0.75735636... 
max nnn 

This curve has an energy which is only 93.42% of that of the semicircle. We have found a curve which has a 
smaller value of 8 than our two-arc solution. 

Can we do better still? 



Multi-arc Approximation 

Consider a smooth curve constructed out of n circular arcs (see Figure 4). Let the radius of curvature of the 
piece turning through die angle from a- to a. +1 be r.. We note that die total arc length, if ', and the integral 
of the square of the curvature, 6 ', are given by 

n-l n-l 

*'= 2 (« i+ l-«l> r l = r n-l a n+2 <'i-i- r 1> a 1' 

i-0 1 n-l i-1 L i i-1 

where a Q = and a n = tt/2. We also have to compute the width, 1f , and height, % ', 

n-l n-l 

W=5J ^(cosa. -cos a 1+l ) = r Q cos a - 2 (>Vi ~ r i )cosa i' 

n-l n-l 

D6'=2 ^(sina^-sina^ Vi sina n + 2 (r._ a - r.)sin a.. 
i-0 i s l 

1 o solve our original problem we need to scale whatever curve we obtain so that its overall width equals 2 
instead of 2Hf '. The scaled values are as follows: 

■ g = 8'ir, 

Note that the integral of the square of the curvature is decreased when we make the curve larger without 
changing its shape. 
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Optimuin Multi-arc Approximation 

Our task now is clear: We have to minimize 8 ' 1f by suitable choices of the parameters r. (for /=0,1 ... n -1) 
and a. (for /= 1,2 ... n -1). Actually, we can pick one of the radii arbitrarily, r Q for example, since the whole 
curve will simply be scaled accordingly. The minimization looks difficult at first when one considers the 
complexity of the product 8 ' 1f and its derivatives. It appears necessary to resort to numerical techniques to 
solve for the 2(n -1) parameters. 



Fortunately this is not the case, for if 



^<6'in = o, 



then, by the rule for the differentiation of a product, 

for arbitrary x 0- e. r. and c^). Since the right-hand side is independent of x, it must equal a (positive) 
constant, c 2 say. Thus we find that 

/-\ ~"37"/"aT =:c ^ for / = 0,1 ... n -1. 



We need the following derivatives now, 



98' 
3 



f: = -[r-7 i :] fo,,/ = 1 - 2 - B - 1 ' 



.1 1-1 



3 g ' K+i" a i> . A1 
_ = _ for r =0,1... /i -1, 

B^"^ ( r i-i ~ r i) sina i for/= 1,2... /i -1, 

1 

g^- = (cos a. - cos o i+1 ) for / = 0,1 ... n -1. 



Using these derivatives in the equations above we obtain 

r. r 1M sin a. = c 2 for / = 1,2 ... /i -1, 



cosa i+1 - cos a. 

a i + l a i 
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It is easy to verify that in the case that n-% we obtain the same equations as before, provided we introduce 
the additional constraint V = 1 or 

r -(r - : r 1 ) cos 04 = 1. 

Note that, if c is known, a simple procedure will give us all of the parameters. Let r Q = 1, say, then the second 
equation can be used to find a y ( This non-linear equation has to be solved numerically.) The first equation 
then allows one to solve for r y Knowing r v the second equation allows one to find a v and so on. If the 
value of c is correct, the process will terminate with a n -7r/2, The correct solution can be found by 
searching for the appropriate value of c. This is very much simpler than a direct search on the 2(n -1) 
parameters. 



Some Helpful Relationships 

A number of interesting observations can be made now about the multi-arc solution. First of all, the "energy v 
66 ' in an individual arc is directly proportional to the projection 8 V of this arc on the x-axis, since 

(a. x1 - a.) 

i 
and 8W = r. (cosa i - cosa i+1 ). 

So we have SS 78 HT = 1/c 2 , 

and we already know, of course, that 6 7 V = 1/c 2 . 

Next, notice that the projection, 8% ', of the (/' + l)-th arc on the y axis equals 

fi06' = r i (sina i+l -sina 1 ) = c z r^--- L -]. 

i+l i-1 

The height at the tip of the /- th arc then is 

% .' = c 2 F 1 + ™ - -^1 for / = 1,2 ... n -1, 
Lr. r^ r Q J 



since 



% 1 ' = r sina t = c z /r r 



'n-1 



Three, Four and Five Arcs 

The optimum solution for three arcs gives 
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% = 1.399926... 
6 = 1.456879... s; 0'.9274780*(w/2) 
S= 1.929128... S= 1.228121*(w/2). 
and shows us that the ellipse is not optimal after all. For four arcs we find 

% = 1.462089... 
8 = 1.448212... SS 0.9219604*(tt/2) 
y = 1.987501... X 1.265282*(w/2), 
and for five arcs we get 

% = 1.500993... 

g = 1.443930... S 0.9192345*(ir/2) 

y = 2.024437... 2Z 1.288796*(vr/2). 

These solutions are shown in Figure 5. We see that 8 is dropping more and more slowly, while f is growing, 
as is %. 

For five arcs, the parameters for the unsealed curve (r Q = 1), are as follows 

r x = 0.500258... ^ = 0.078707... 

r 2 = 0.334497... a 2 = 0.237276... 

r 3 = 0.253163... a 3 = 0.483045... 

r 4 = 0.207337... a 4 = 0.847067... 

r 5 = 0.189705... a g = 1.57079... 

Perhaps we can guess the true minimum energy curve from the numerical data obtained so far. The 
parameters seem to roughly fit into a pattern like 

1 0+1) i 2 n(n+l) 

In this case ( a a ) = w _J±I_ 

i+i t' n(n+l)' 

so the arc lengths are 8? ' = r.(a .^, - a .) = , ff | ,, . 

i i+i v n(n+\) 

/"""N That is, the arcs all have the same length, and curvature increases linearly along the curve. 
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Now wc could sum the scries for arc length and energy (easy) and for height and width (hard). Then we 
would discover that 8 = 8 ' W decreases with «, and we could find its limit as the number of arcs tends to 
infinity. Instead, we proceed directly to the curve obtained in the limiting process. 



The Cornu Spiral 

The curve which has curvature varying linearly with arc-length is called the Cornu Spiral (or Ruler's Spiral) 
[pg. 190, 2; pg. 190, 14]. It can be defined using the two Fresnel integrals [pg. 820, 10; pg. 930, 11; pg. 300, 
12]. 

C(s) = f $ cos (-fydt, 

S(s) = f $ sin (£( 2 )dL 

If we let x= C(s) and y = S(s\ 

we obtain a curve starting at the origin and curling upwards in the first quadrant. We note that 

x = cos (ys 2 ) and y = sin (^X 

x = - ms sin (ys 2 ) and y = 7rscos(ys 2 ). 

This verifies that s is the arc-length along this curve, since 

jc 2 + y 1 = 1, 

and that the curvature varies linearly with arc length, since 

y x - 3c y = *n s. 

The part of the spiral of interest to us here extends to the right up to the point where the curve becomes 
vertical, that is, x- 0. This is the point where s~ 1 and 

x = CU) = 0.7798934... and y = S(l) = 0.4382591... 

The energy in this portion of the spiral is just 

8'= [ X K z ds= f 1 (vsfds = ?r 2 /3. 

We now build a smooth curve connecting the two points (-1,0) and (+1,0) by scaling, rotating, and shifting 
this tendril as shown in Figure 6. In die right hand quadrant we use 

x - 1 S{1) and y - m 

(The rest of the curve is obtained by reflection about the yaxis.) 
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We then find that 

Dfi = C(l)ASXl)= 1.779525... 

6 = (w 2 /3)6Xl) = 1.441814... S0.9178877*(w/2), 

S = 1/5(1) = 2.281755... « 1.452610*(tt/2). 

The curve constmcted out of a portion of the Cornu Spiral only has 91.78% of the energy of the semicircle 
and is thus the best curve so far. 

But, can we do better still? 

Six Arcs and More 

Unfortunately, the Cornu Spiral is not optimal either, as one sees by considering the best six-arc solution for 
which 



/^ 



For eight arcs, 
for sixteen arcs* 
for thirty-two arcs, 
and for sixty-four 



6 = 1.441508... S 0.9176931*(tt/2). 

6 8 S 0.916097*(t7/2), 
S 16 ~ 0.914532*(tt/2), 
e 32 ~ 0.914285*(t7/2), 
6 64 ZZ 0.913953*(tt/2). 



These solutions are shown in Figure 7. 

It seems that the total energy is approaching some limit, near 91.39% of that in the semicircle. Some of these 
results are summarized in Tables I and II. 

TABLE I 



8/(97/2) 



*/(ir/2) 



% 



1 


1.0 


1.0 


1.0 


1.0 


2 


.9415383 


1.160608 


1.284916 


.9245284 


3 


.9274780 


1.228121 


1.399925 


.8950140 


4 


.9219604 


1.265282 


1.462089 


.8794932 


5 


.9192345 


1.288796 


1.500993 


.8700129 


6 


.9176896 


1.305012 


1.527622 


.8636634 


7 


.9167300 


1.316870 


1.546987 


.8591353 


8 


.9160932 


1.325918 


1.561701 


.8557555 



^■-y 
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TABLE 11 




n 


8/(w/2) 


*/(ff/2) 


% 


min 


16 


.9144692 


1.358673 


1.614492 


.8442353 


32 


.9140406 


1.375704 


1.641625 


.8389106 


64 


.9139305 


1.384403 


1.655392 


.836516 


128 


.9139025 


1.388789 


1.662308 


.835461 


256 


.9138955 


1.390998 


1.665786 


.834998 


512 


.9138938 


1.392110 


1.667534 


.83473 


1024 


.9138934 


1.392671 


1.668417 


.83467 


2048 


.9138932 


1.392984 


1.668908 


.83463 



We can get better and better approximations to the optimum curve, provided we also carry out computations 
with more and more significant figures. Note, by the way, that while £ varies little once n is reasonably large, 
y and % continue to show appreciable changes. This is a reflection of the fact that some distortions of the 
optimum curve produce only small changes in the total energy. 



Some Observations About the Optimum Curve 

The multi-arc approximation tends to the optimum curve in the limit as n tends to infinity. So we can learn 
some properties of the optimum curve from what we have so far. First of all, 



from 



we get 



r.r J ._ 1 sin a^ = c 2 



C 2 K = COS *|>, 



where k is the curvature and ^ is the angle which the curve makes with the x-axis. The constant c only affects 
the size of the curve, not its shape. We cannot determine it at this point. 



From 



we get 



cosa m _cosa. 2 

- r. — ~ — = c , 

.sin [(a.. --a J/2] x ' ? 



which in the limit again leads to 
We also obtain 



C 2 K 2 = COS $. 

dx 



from 



— - Mr 7 - 

6W ~ Uc ■ 



Now 



g = JV<fr= Jk 2 -/ 1 + (dy/dxf dx, 
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^ 



r^s 



so 



K 2 -y/ 1 + (dy/dxf = 1/c 2 . 



Further, since, -^ = tan J/, 

we again obtain, c 2 * 2 = cos $. 

Each of these approaches leads us to the same simple differential equation for the curve. 

Finally, from %' = C 2r± + _L-_±1 

Lr i r i-l r J 

we get in the limit - k = >>/2c 2 . 

That is, the curvature varies linearly along the axis of symmetry of the optimal curve. Substituting for k we 
also derive 

cos \p = (y/2c) 2 . 

Note that since the optimal curve bends downwards, its second derivative is negative. This is why, by the 
usual sign conventions, curvature too is negative. Thus we will use the equation 

— C K = V COS\p 

between 4> = + 77/2 at the left end and \p - - tt/2 at the right end of the curve. 

Differential Equation for the Curve 

The curvature is the rate of turning as one goes along the curve, that is, 

K ~ ds> 

so - c-% = V cos $, 

and consequently, s = - c j / r . 

The substitution cos \p = £ leads to a denominator of the form 

V*(l-*)(l + «), 

while the substitution cos \p = £ 2 leads to a denominator of the form 

VTT7. 

/ ^ sss \ ) In each case we are dealing with the integral of a rational function of £ and the square root of a cubic or 

quartic polynomial in £, This means that the answer can be expressed as an elliptic integral [pg. 16, 9; pg. 833, 
10; pg. 904, 11; pg. 589, 12]. 
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Wc can actually just look up the result directly [using pg. 154, 1 1 ; also see Appendix] and find the solution 

^,/2cF(cos^V r coT^, 1/72), 

where F is the incomplete elliptic integral of the first kind. 1 The constant of integration has been chosen so 
that f=0 corresponds to s=-0. Hie result came out positive because the integration goes from i/< = to 
negative values of */>. For the half of the optimal curve in the negative quadrant, a minus sign must be 
attached. Some readers may notice the similarity of the solution found above, to the equation for a pendulum 
[pg. 28, 9], swinging from -mil to -\-mll (where i> is the angle from the vertical, while s corresponds to the 
time). 

Wc now have the equation of the curve sought after in Whcwell form [pg. 4, 14], namely arc length as a 
function of tangential angle. We can immediately also rewrite it in Ccsaro form [pg.4, 14], namely as a 
relation between arc length and curvature: 

s = VlcFicos'H-CKll/SD. 

Both of the forms given above are intrinsic equations for die curve [pg. 40, 2]. 

We can easily compute the length of the curve from an initial point at the top, where $ = 0, to the point on the 
jc-axis, where $ = - tt/2, since 

/< (tt/2, 1/71) - K(1/V1) 

where A' is die complete elliptic integral of the first kind. Now [pg. 909, 11], 

m/ST) = *(sin(T/4)) = v^ 1 j *_ ^ = gffig, 

where T is the gamma function [pg. 821, 10; pg. 933, 11; pg. 255, 12]. There is an infinite product [pg. 938, 11] 
forr(l/4) 4 



r(l/4) -16* U (4yt+1) 2 [( ^_ 1)2 _ 1] 



which gives us the numerical value 



F(l/4) = 3.6256099082.. 
cr(l/4) 2 



The arc length is finally, J = ^ 7 ^ • 



1. Note that some authors [pg. 833, 10] list the arguments of the incomplete elliptic functions in the reverse 
order of that shown here. 
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Cartesian Form of the Solution 

For many purposes it is more convenient to express the solution as a relationship between the x and y 
coordinates. 

We note that —■ = cos ^ and J = sin */>, 

and remember that - c~r = V cos \j,. 

By the chain-rule for differentiation 4f = 4* I ^ = - c/comF, 

dy ds* ds T ' 

and ^ = g / ^ = _ csin + / ^« 7 . 

The latter equation is easy to integrate using the substitution z = cos i£. 

dz 



y 



= c ff- 2 = 2cV ~ z > 



so that (y/2c) = V cos \J>, 

where the constant of integration was chosen so that ^=0 when ^= - tt/2. The reader may also recall that in 
f\ the previous section this result was found directly as the limit of height of the multi-arc approximation. 



The integral for x is a little harder, x = c j Vcosi^ #. 



This can be expressed as the difference of two incomplete elliptic integrals [using pg. 156, 11; also see 
Appendix], 

x= v^c[2i^cos" 1 Vlos7,l/v r I)- Fico^Vwsy.l/VDl 
where the constant of integration was chosen so that x=0 when i// = 0. Finally then, 

x = VI c[2 E (cosine), 1//1) - F(cos~ V2c), 1//!) J. 
An alternate way to obtain the same result is to note that x and y are related by the differential equation 

dv _ V 1 - iy^cf 

dx (yllcf 



1. Note that E denotes two things in this paper. When the letter appears with two argument it signifies the 
incomplete elliptic integral of the second kind. 
^n 2. Note again that some authors [pg. 833, 10] list the arguments of the incomplete elliptic functions in the 

reverse order of that shown here. 
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since X = tan *// = 



"•* COS ip 



In any case, for j = 0, 



x=/2c[2m//2)-«l//2)], 
where K and £ are the complete elliptic integrals of the first and second kind respectively. 
Using Legendre's identity [pg. 25, 9; pg. 836, 10; pg. 90J11; pg. 591f 12] 

E(k)K(k') + E(k') K(k) - K{k)K{k') - w/2, 
where k z + k' 2 = 1, 

we get [pg. 25, 9] 2 £(1/72) = [ ir/2 + # 2 (1/ VI) ) / K(llVl) , 

so that [ 2 £ll//2) - K(ll VI) ) = <w/2) / AW VI) . 

(2tt) 3/2 
So the width of the curve is W = c p/yA2 » 

and, if we want H = 1, we must have 

c = ;^{^2 = 0.8346268416... 
(2w) 

The height of die curve then comes to 

36 = 2c = 1.669253683... 

(compare to the two-arc approximation). The minimum radius of curvature, the inverse of die maximum 
curvature, is 

r . = ] 7 = c = 0.8346268416... 

TTius a circle tangent to the curve at the top is also tangent to the x-axis. 

1 F(l/4) 4 
The arc length comes to if = 5 1 \jT > 

or J = 2.188439615... « 1.393203929*(w/2). 

Finally, from ^r = 1/c , 

we get 8 = 1/c 2 , 



and so 6 = ^1_ = 1.435540022... « 0.9138931623*(w/2). 



Note that 6 $ = w. The curve of least energy is shown in Figure 8. 
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Extension of the Curve 

So far we have considered a finite segment of the optimum curve, extending from a point of zero curvature 
(where i// = -tt/2) through a point of maximum curvature (where i// = 0) to a second point of zero curvature 
(where ^ = -f it/2). Can the curve be extended beyond these points? 

It is clear that ^ must remain in the range [—tt/2, +tt/2] so that the square root of its cosine remains real. 
To continue the curve then, the sign of the curvature must change; we must choose the other sign for the 
square-root. The new segment we obtain has the same shape, of course, as the segment we have found 
already, just inverted. 

The segment we have used so far is just a piece of the infinite periodic wave shown in Figure 9. In Figure 10 
we see several curves which correspond to stationary values of the integral and pass through the specified 
points with the desired orientation. The one on the left is the one which corresponds to a global minimum of 
the energy. The curves containing n half-cycles have an energy n 1 as large as the one containing a single 
half-cycle. 

The curve of least energy passing through two given points with specified orientation is just a portion of the 
general curve, suitably translated, rotated, and scaled. This is illustrated in Figure 11. The rotated, translated 
and scaled curves form a four parameter family. 



Variational Approach 

^n We are trying to find the curve for which 8 = f k 2 ds 

is minimal. This integral can also be written in the form 

g = JV[1+00 2 ] 1/2 rfx, 



or, since K = 



JL 



[1 + iyff n ' 
til 



This is of the form f<jf ( x , y, y\ y") dx, 

and the calculus of variation [pg. 119, 15; pg. 190, 16; pg. 198, 17] teaches us that for a stationary value of the 
integral, 



Act ,.jLL 



5v-^v' + l2 ? v" = 0, 



where f y , 9y, and *F y » are the partial derivatives of % with respect to y, y, and /' respectively: 



r> 
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V 


0, 

5y(y") 2 
' [l + (y'f) 1 ' 2 


V ! 


2y" 


" [i + 00 2 ] 5/2 " 


-IK' 


-lvH< 



Since 9 = 0, 

and integrating, we get -9 y ' + ^9y< = A 

where J is an arbitrary constant. In the above we have closely followed the approach taken by Mehlum 
[pgs. 157&189, 2; pg. 43, 



A<x 2y'" 10 y'(y') 1 

</* y" _ [l + (y'ff n [l + Cvf] 7 



AlSO —^v" = ri , / <\2i5/2 ~ ri _l f./\2l7/2 < 



2/" 5y'(/) 2 . 

so we get [x + ^ 2]6/2 " u + (Jyjy/z - >*■ 



Now 



dy" dy" £ dV ,,, , ,i 



so >- =y -%r- 

j j- • ^ 1 / "\2 " dy" 

In addition dy'2^' ~ y dy ' 

A d 1 5y 

and ^[l + ^) 2 ] 5/2 ~~[l + (y') 2 ] 7/2 



Using these results in the equation above, 



A^Jzl 



dy\l + (y?? n 

iy'f 

and, integrating, we get n . < ^2i5/2 = Ay + B 9 

where /? is a second arbitrary constant. 

Returning for a moment to the integral of the square of curvature, we see that 
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fn z ds = J [Ay + B] dx = J Ady+ f B dx 
for a curve which starts at O , y ) and ends at (x^ , j^). 

now |1+ ( $ r = «yi +</)'. 

so k 2 yHTfTyf = >< y + a 

Also, <k = & f <k = p/i 

dx ds 1 ds y ' 

SO K 2 -^ 1 + OVi) 2 = yj (y/i) + #, 

or /c 2 = ^. -f M, 

since i 2 + j^ 2 = l. 

Nowif tan* = ^, 



then i = cos i/> and j> = sin \p 

Remembering that 






we finally see that d = ± ^ A sin ^ + 5cos ^ • 

Letting yj = 1/c 2 cos cp and Z? = - 1/c 2 sin <p 

weget f 7fc = ± ^ sin (^ - ?)• 

The scale of the curve is dependent on the parameter c, while the rotation in the xy-planc is dependent on the 
parameter <p. 

Altogether we have a four parameter family of curves, since we can also choose an initial point and a direction 
for the curve. Conversely, we can find a single curve out of this family which passes through any two points 
with orientation specified at both points. 

By the way, if we let the line from the initial point (x Q9 y ) to the final point (x^) have length r and 
direction 6, then 



JV& = (r/c 2 ) sin (0 -<p). 
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Summary 

Wc have found the simple equation 

± ck - V cos (lp - (p) , 
and solved it to find the equation of die curve of least energy in Whewell form 

5= ±v^cF(cos" a /cosTF^,l/V^). 



We also developed a differential equation for the curve in Cartesian coordinates aligned with the axis of 
symmetry, 



V 1 - iyncf 



dx (y/2cf 

The solution of this equation, for given initial conditions, led to 

x = /2c [ 2 E(cos^(y/2cl 1//2) - F(cos" V2d 1//5)] . 

Wc considered the curve of least energy connecting the point (-1,0) to the point ( + 1,0) with vertical initial 
and final orientations. This curve has minimum radius of curvature 

F(l/4) 2 
C ~(2tt) 3/2 ' 



rises to a height 3G = 2c, 



4 



i A t i ra/4r 

has arc length 2 o ^ 2 ' 

e (2tt) 3 

and energy 8 - r(1/4) 4 > 

or about 91.39% of that of the simple semicircle approximation. 

We have also given a method for finding approximations, consisting of circular arcs, to the curve of least 
energy. 

Note that the curve found here is extensible [7] in the sense that if the least energy curve with orientation a at 
A and orientation at B passes through the point C with orientation y, then the segments from A to C and 
from C to B are themselves least energy curves. (This is not true of the two-arc approximation [8]). As a 
result, such a curve can be computed by a simple, locally connected network. 

We have not shown how to find the particular member of the four parameter family of curves which passes 
through a specified pair of points with specified orientation. Presumably determining the axis of symmetry 
of the curve would be a helpful first step in this direction. In practical applications the multi-arc 
approximation method may be suitable in dealing with this problem. We have not discussed how one might 
compute the curve of least energy passing through three or more points. Here there is no constraint on the 
direction, but the curvature must be continuous. Nor have we touched upon the extension to curves and 
surfaces in three dimensions, a topic which Mehlum addresses [pg. 62, 6\. 
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Appendix 

In this appendix we determine two integrals needed in the body of the paper. 



y^^s 



To evaluate the integral L = f -7^— ? 

1 Jo Vcos^' 

we substitute £ 2 = cos i/>, 

V cos^ d£ 



and obtain / = — 2 f 

1 *>i 



Next, substituting t 2 = 1 - £ 2 , 

r r* r V 1 - cos ^ <# 

weget / = /2 I r , . 

1 J ° Vl-/ 2 7l-/ 2 /2 

and since sin" * V 1 - cos ^ = cos" a V cos t/>, 

we finally get J^ ygL^ = /2 F^os" 1 /S^, 1//5). 



f*\ 



To evaluate the integral l z = f Vcos^ rf«|/, 

we again substitute £ 2 = cos ^ 

, i . r r V cos i£ £ 2 

and obtain /, = -2 — , <ff 

m i 2 i + g 2 - 1 VT7I 

NOW / = 1 > ■ = Y , - ; . 

The integral thus can be split into two parts, the second of which we have already evaluated. 



2 



-22- 



To evaluate L - L = - I / = a$ , 

2 l Jl V 1 - e 



2 _ i t 2 



we substitute as before / = 1 - £ 



- rVl -cos* V ] - / 2 /2 

and obtain / -/ 1 =2V2I ^^====^ dt , 

2 1 j y j , 2 



f 



andso 7 2 - ^ = 2-/1 /< (cos" 1 /"^?, 1//2) , 



since 



sin^V 1 - cos * = cos" 1 -/ cos* 



Finally then, 



f /cos^# = V r 2[2/;(cos*' 1 v r coF?,l//2)- F(cos" a v r cos?, 1//2)]. 



The corresponding definite integrals from * = to $-<n/2 can be expressed in terms of complete elliptic 
integrals of the first and second kind [pg. 91, 9]. 
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Figure 1] : Optimum curve passing through two given points with given orientation. It is a member of a four 
parameter family of curves obtained by translating, rotating, and scaling the particular curve 
passing through the points (- 1, 0) and ( + 1, 0) with vertical orientation. 



